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Abstract 

We present an analysis of mock spectral observation of warm-hot intergalactic 
medium (WHIM) using a constrained simulation of the local universe. The simulated 
map of oxygen emission lines from local WHIM reproduces well the observed struc- 
tures traced by galaxies in the real local universe. We further attempt to perform 
mock observations of outer parts of simulated Coma cluster and A3627 adopting the 
expected performance of BIOS (Diffuse Intergalactic Oxygen Surveyor), which is pro- 
posed as a dedicated soft X-ray mission to search for cosmic missing baryons. We 
find that WHIMs surrounding nearby clusters are detectable with a typical exposure 
time of a day, and thus constitute realistic and promising targets for BIOS. We also 
find that an X-ray emitting clump in front of Coma cluster, recently reported in the 
XMM-Newton observation, has a counterpart in the simulated local universe, and its 
observed spectrum can be well reproduced in the simulated local universe if the gas 
temperature is set to the observationally estimated value. 
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1. Introduction 

It is widely accepted that our universe is dominated by dark components; 23 percent in 
dark matter, and 73 percent in dark energy (Spergel et al. 2003). More surprisingly, even the 
remaining 4 percent, cosmic baryons, has largely evaded the direct detection so far, i.e., most 
of the baryons is indeed dark. While those cosmic missing baryons may consist of compact 
stellar objects (white dwarfs, neutron stars and black holes), brown dwarfs, and/or diffuse gas, 
recent numerical simulations indicate that they largely take a form of the warm-hot intergalac- 
tic medium (WHIM) with temperature of lO^K < T < lO'^K (e.g., Cen & Ostriker 1999a). In 
fact cosmological hydrodynamic simulations consistently point to the fact that WHIM traces 
the large-scale filamentary structure of dark matter distribution more faithfully than hot intr- 
acluster gas (T > 10''K) and galaxies both of which preferentially resides in clusters that form 
around the knot-like intersections of the filamentary regions. This implies that WHIM carries 
important cosmological information in a complementary fashion to distribution of galaxies (in 
optical) and of clusters (in X-ray). 

Unfortunately the X-ray emission of WHIM via the thermal bremsstrahlung is very weak, 
and its detection has been proposed only either through Ovil and Ovill absorption features in 
the QSO spectra or the possible contribution to the cosmic X-ray background in the soft band. 
While there are several observational reports to have detected signatures of WHIM via oxygen 
absorption lines (Nicastro et al. 2002; Fang et al. 2002; Mathur et al. 2003) and via excess 
X-ray emissions (Kaastra et al. 2003; Finoguenov, Briel, & Henry 2003) among others, they are 
not well suited for unbiased exploration of WHIM that is important for cosmological studies. 
In order to identify the existence, and explore the nature, of elusive cosmic missing baryons 
in a more unambiguous manner, we propose a dedicated soft X-ray mission, DIOS (Diffuse 
Intergalactic Oxygen Surveyor; see Ohashi et al. 2004) which will carry out a direct and 
homogeneous survey of WHIM. In a previous paper (Yoshikawa et al. 2003; Paper I hereafter), 
we examined in detail the detectability of WHIM on the basis of the cosmological hydrodynamic 
simulations. We indeed showed that the unprecedented energy resolution (~ 2eV) of the X-ray 
Spectrometer Array (XSA) on-board DIOS enables us to identify WHIM with gas temperature 
T = 10^~^K and overdensity 6 = 10 — 100 located at redshift z < 0.3 through emission lines of 
Ovil and Ovill. 

The purpose of the present paper is to consider the detectability of WHIM surrounding 
nearby clusters {z < 0.03) using the constrained numerical simulation to reproduce the local 
universe (Dolag et al. 2003). This approach is useful in setting the observational strategy of 
the targeted observations around the previously known structures. This is a complementary 
and equally important observing mode to the blind survey that is implicitly kept in mind in 
Paper I. 

It should be noted here that after submitting the first manuscript of the present paper, 
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we were informed of an earlier paper discussing exactly the same topic by Kravtsov, Klypin 
& Hoffman (2002). They used the constrained simulation of the nearby universe based on 
(Klypin et al. 2003) the MARK III survey of peculiar velocities, and addressed the observational 
signatures including emission by Ovil and Ovill. Since they do not specifically bear in mind the 
high-resolution spectroscopy which will be feasible with BIOS, they do not study the detailed 
emission spectra that we will present below. Nevertheless the overall qualitative features that 
they found are in good agreement with our results. 

The rest of the paper is organized as follows; section 2 describes the local universe 
simulation data that we use for mock observations. In §3, we summarize how to compute the 
soft X-ray spectrum of the local universe with the expected performances of BIOS following 
Paper I. Several examples of mock WHIM observation are shown in §4. We focus on the 
mock observation around the outskirts of the simulated Coma cluster in §5, and show that our 
simulated spectrum reasonably reproduces a recent observational report of the XMM-Newton 
discovery of an X-ray clump in Coma by Finoguenov, Briel, & Henry (2003). Finally, section 
6 is devoted to the conclusions of the paper. 

2. Simulation of the Local Universe 

Throughout this paper we analyze a dataset generated by a cosmological hydrodynamic 
simulation which attempted to reproduce our local universe. The simulation is based on the 
previous work by Mathis et al. (2002) which intentionally constructed the initial conditions 
making use of the IRAS 1.2-Jy survey (Fisher et al. 1994; Fisher et al. 1995) out to cz < 
12000km s~^. They assumed a ACDM (A-dominated cold dark matter) model with the density 
parameter = 0.3, the dimensionless cosmological constant = 0.7, the Hubble constant 
h = 0.7 (in units of 100 km/s/Mpc), and the rms density fluctuation ag = 0.9 (top-hat smoothed 
over a scale of 8/i~^Mpc). First the observed galaxy density field was Gaussian smoothed on 
a scale of 5/i~^Mpc. Then it was evolved (quasi-)linearly back to the initial redshift of the 
simulation, z = 50. Its overall amplitude is normalized so that the resulting rms mass variance 
matches the adopted value of as- This was used as a Gaussian constraint for an otherwise 
random realization of the ACDM cosmology. Centered on the Milky Way the size of simulation 
box is 240/i~^Mpc, with a high resolution region covering a sphere with 80/i~^Mpc of radius 
which corresponds to the region constrained by the observations. Evolved forward in time, 
Mathis et al. (2002) demonstrated that this simulation very well reproduces the local universe 
today in several respects; applying semi-analytic galaxy formation scheme to their dark matter 
only simulation very well reproduces the local galaxy population. Density and velocity fields 
obtained from synthetic mock galaxy catalogues reproduce very well the observed counterparts. 
Also some of the most prominent halos in the simulations may be identified and indeed agree well 
in positions and masses with the observed galaxy clusters. Further details of the construction 
of the initial condition are found in Mathis et al. (2002); see also independent earlier work by 
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Klypin et al. (2003) and Kravtsov, Klypin & Hoffman (2002). 

Dolag et al. (2003) repeated the simulation including gas physics with smoothed par- 
ticle hydrodynamics (SPH) technique. They extended these initial conditions by splitting the 
original high resolution dark matter particles into gas and dark matter particles with masses 
of 0.48 X 10^h~^MQ and 3.1 x lO^/i'^M©, respectively. The gravitational force resolution of the 
simulation was set to be 7.5/i~^kpc which is also the resulting mean inter-particle separation 
of the SPH particles in the dense centers of the simulated galaxy cluster. The most massive 
cluster in this simulation is thereby resolved by nearly one million particles. 

This simulation was carried out with GADGET-2, a new version of the parallel TreeSPH 
simulation code GADGET (Springel, Yoshida & White 2001). It uses an entropy-conserving 
formulation of SPH (Springel & Hernquist 2002), and was supplemented with the treatment 
of magnetic field in the framework of ideal magneto-hydrodynamics as described in Dolag et 
al. (2002). Since the primary purpose of the simulation was to study the propagation of the 
ultra-high energy cosmic rays in the extra-galactic magnetic fields, the simulation did not in- 
clude the effects of radiative cooling of gas, galaxy formation, energy feedback and metal enrich- 
ment. Nevertheless the simulation data provide useful mock samples for targeted observations 
of WHIM; for example, according to numerical simulations presented in Dave et al. (2001), 
the effect of gas cooling can change the dumpiness of WHIM. The difference between the re- 
sults with and without radiative cooling, however, is within a factor of two. Therefore, it may 
change the results at some places in a quantitative manner, but as far as the detectability 
and the overall distribution of WHIM (smoothed over the angular resolution of BIOS, for in- 
stance) are concerned, such details are not supposed to be essential. On the other hand, the 
detectability of the WHIM via oxygen emission lines is sensitively dependent on the assumed 
metallicity model. So we adopt an empirical model (eq. [1] below), which is shown to repro- 
duce the observed typical metallicity of intracluster medium (Cen & Ostriker 1999b; Aguirre 
et al. 2001). 

For our analysis, in order to avoid a possible unrealistic effect due to the numerical 
boundary, we restrict ourselves to a region within 75/i~^Mpc from our Galaxy. Figure 1 shows 
the Hammer-Aitoff map of baryonic matter distribution in the supergalactic coordinate. As al- 
ready shown in Mathis et al. (2002) and Dolag et al. (2003), the simulated distribution of matter 
reproduces well many known nearby structures, including the Great Attractor region, the Great 
Wall structure, the Pisces-Perseus supercluster, and the Local Void, as well as many known 
rich clusters such as Virgo, A3627 and Coma. Additionally, the temperature and mass of the 
identified halos in the simulation roughly agree with the observed values of their counterparts 
in the real universe. Table 1 compares the simulated and observed values of celestial position 
in the supergalactic coordinate, redshift, mass and emission weighted temperature for Coma, 
Virgo, A3627, Hydra, Centaurus and Perseus clusters. The quoted masses in Table 1 denote 
the virial mass within -R200; the radius inside which the mean density is equal to 200 times the 



4 



critical mass density for Coma, Virgo, Hydra, and Perseus, while we consider mass within 40' 
for A3627 and 50' for Centaurus cluster. They are consistent with the observed values, implying 
that the simulation data are reliable enough for our present purpose. Nevertheless we would 
like to emphasize that despite the very promising comparison on large scales, many details in 
the simulations are not fully constrained by the construction of the initial conditions. That is 
specially true for the filamentary structure, which is certainly correct in representing their pres- 
ence in a statistical manner within the supercluster structures, but their exact position within 
these complexes cannot be compared with the real observations. Therefore predicted emission 
from these structures within the superclusters can be well compared with observations, but one 
should expect that the exact location within these superclusters will not be reflected in real 
observations. 

Within the supergalactic coordinates, we have many galaxy clusters and superclusters 
near the supergalactic equatorial plane. The upper left and right panels of Figure 2 show the 
density and temperature field of a 9/i~^Mpc thick slice parallel to the Super-Galactic plane, 
extending from SGZ = -4.5/i-^Mpc to SGZ = 4.5/i-^Mpc. 

3. Oxygen Line Emission from WHIM in the Simulated Local Universe 

The goal of this paper is to examine the detectability of WHIM in the local universe 
through its oxygen line emission. Naturally the expected flux sensitively depends on the metal- 
licity of the WHIM or more generally of the intergalactic medium (IGM), which is very uncertain 
both theoretically and observationally. In addition the simulation data that we use do not take 
account of the effect of galaxy formation and its succeeding energy feedback and metal enrich- 
ment. Thus we have to assume a phenomenological metallicity model in subsequent analyses. 
In this paper, we adopt the following density- dependent model: 

Z/Zo = Min [0.2, 0.02(pgas/pb)°-'] , (1) 

where Zq denotes the solar abundance. This is motivated by the radiation pressure ejection 
model of Aguirre et al. (2001) where the stellar light exerts radiation pressure on the interstel- 
lar dust grains and expels them into the hosting galactic halos and the ambient IGM. This is 
the same as Model IV in Paper I except that we put an upper limit 0.2Zq so as to avoid the 
unacceptably large metallicity at cluster center (Gen & Ostriker 1999b). In practice, however, 
the upper limit does not affect the result for WHIM that is of our primary interest. The lower 
left panel of Figure 2 indicates the resulting metallicity map on the Super-Galactic plane, illus- 
trating that the prominent filamentary structures and supercluster regions exhibit metallicity 
of logio(Z/Zo) > -1.5. 

The proper identification of the oxygen emission lines from WHIM is technically very 
challenging with the existing X-ray missions, and requires a small but dedicated soft X-ray 
mission (Paper I). Throughout the paper, we assume the currently proposed specification of 



DIOS (Ohashi et al. 2004) which is summarized in Table 2; high energy resolution of Ai? = 2eV 
around E ~ 600eV, and a nested 4-stage X-ray telescope with an effective area Scs = 100 cm^ 
and field-of-view fipov = Ideg^ (16 x 16 pixels in total). Table 3 shows the detection limit of 
the line emission at photon energy E = 600 eV for DIOS in the cases of signal-to-noise ratio 
S/N = 10 and 5 assuming three different exposure times, T^^p = 10^, 10^, and 10^ sec. Those 
detection limits are computed from equation (4) of Paper I assuming that the flux of the cosmic 
X-ray background (CXB) is 3 x 10~^ erg s~^ cm~^ sr~^ keV~^ between 0.5 and 0.7 keV. In what 
follows, we adopt the nominal detection limit of 1 x 10~^^erg s~^ cm~^ sr^^ both for Ovil and 
Ovill lines. 

As in the Paper I we calculate the Ovil and Ovill line emission assuming that WHIM 
is under the coUisional ionization equilibrium. In calculating surface brightness of the line 
emission, we divide the whole simulation volume into cells, each of which has colatitudal and 
azimuthal extension of ASGL = 0.5° and ASGB = 0.5° on the celestial plane, respectively, 
and has an equally spaced redshift interval (75 bins from z = to z = 0.03). The line surface 
brightness of Ovil and Ovill for each cell is given by 

AAi = [d„(zt)]^A:j, (3) 

where the summation is carried out over all gas particles (labelled i) within the cell. In the 
above expression, X is the hydrogen mass fraction (we adopt 0.755), nip is the proton mass, 
e(T, Z) is the oxygen line emissivities for gas temperature T and metallicity Z in units of the 
power input normalized to the electron density as defined by Mewe, Gronenschild, & van den 
Oord (1985), fc{T) is the electron number density fraction relative to the hydrogen, d^iz) is 
the angular diameter distance to redshift z, Au is the solid angle of each cell, and Tj, Zi, pi, and 
rrii denote temperature, redshift, mass density and mass of the i-th gas particle, respectively. 

We note that since we assume that WHIM is in the coUisional ionization equilibrium 
and not affected by photoionization by CXB and UV background radiation, both the electron 
number density fraction /e and the oxygen line emissivities e are independent of gas density. 
According to Kang et al. (2004), more than 95% of the Ovil and Ovill line emissions come from 
WHIM with temperature T > lO^K even under the presence of photoionization background. 
Actually, under the ionization equilibrium, the emissivities e of Ovil and Ovill at T > lO^K 
is little affected by the photo-ionizing background. Thus, as long as ionization equilibrium 
is assumed, the results presented in this paper are not significantly altered even if the effect 
of the photo-ionization background is considered. For WHIM with lower temperature and 
lower density, however, we have to take account of the photo-ionization effect as well as non- 
equilibrium treatment which we hope to report separately. 

Figure 3 shows the simulated all-sky maps of surface brightness for Ovil (574 eV; top) 
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and Ovill (653 eV; middle) line emissions in comparison with the conventional X-ray map (0.5- 
2.0 keV; bottom). At a photon energy of i? = 600eV, the dominant interstellar absorption is due 
to the Galactic oxygen as well as other metals. The oxygen column density is usually assumed to 
be in proportion to the Galactic neutral hydrogen column density. So one can use the observed 
column density of neutral hydrogen A^hi in inferring the degree of the absorption. We consider 
the effect by combining the Galactic Hi map by Dickey & Lockman (1990) and the absorption 
cross section by Morrison & McCammon (1983) where the solar metal abundance is assumed. 
In contrast to the 0.5-2.0 keV X-ray emission, Ovil and Ovill emissions preferentially come 
from small objects like galaxy groups and the outskirts of rich galaxy clusters and superclusters. 
Furthermore, as is clear for simulated Virgo, Ovil (574 eV) emission avoids the very central 
regions of rich galaxy clusters. This is because the emissivity of Ovil emission lines sharply 
drops beyond T > 10^ K (see Fig. 3 of Paper I). Figure 4 depicts the Ovil and Ovill surface 
brightness distribution on the supergalactic equatorial plane (—15° < SGB < 15°) in redshift 
space (in plotting these maps, we use the velocity field of the simulation data, which is smoothed 
over 30 neighbor particles in an SPH manner). A vacant strip extending at SGL = 0° and 
SGL = 180° is due to the Galactic extinction. 

These plots help identify several filamentary structures which have detectable levels 
of oxygen emissions with BIOS; z = 0.01 - 0.015 and SGL = 190° - 200° in front of A3627, 
z = 0.01 — 0.02 and SGL = 280° — 300° in the westward of Pisces-Perseus supercluster. In 
addition, we note an oxygen line emitting clump in front of Coma at redshift z ~ 0.01. We will 
discuss the latter in much detail so as to consider the relevance to the observational report by 
Finoguenov, Briel, & Henry (2003) in section 5. 

If both Ovil and Ovill emissions are simultaneously observed, one can use their line 
ratio as a diagnostics for the nature of the observed WHIM, in particular its temperature. 
Figure 5 shows the contour of correlation between the 0.5-2.0 keV band emission weighted gas 
temperature and the ratio of Ovil (574 eV) and Ovill (653 eV) surface brightness for all the 
regions where both Ovil and Ovill emissions exceed the detection limits (we consider both 
10~^^ and 10~^^erg s~^ cm^^ sr~^ for definiteness) . The surface brightness corresponds to an 
aperture size of ASGL = ASGB = 0.5° {left panel) and 0.25° [right panel), and is projected over 
the redshift interval of Az = 0.0033, which corresponds to spectral energy resolution of DIOS 
[AE = 2 eV) at photon energy E = 600 eV. 

Figure 5 implies that the temperature range of WHIM oxygens detectable by DIOS is 
mainly T = 10^ — 10^ K as already concluded in Paper I. The solid line indicates the relation 
for the collisional ionization equilibrium. Since this is the assumption that we adopted, the 
simulation data should exactly follow the relation if the smoothed cell size is sufficiently small. 
In reality, however, the effective observational resolution is determined by the telescope angu- 
lar resolution and the spectroscopic energy resolution. Thus the inhomogeneous density and 
temperature structure within the cell results in the departure from the assumed relation both 
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statistically and systematically. The systematic deviation from the assumed relation seems 
prominent for lower temperature cells (T < 2 x 10^ K). This is mainly because the oxygen line 
emissions in those cells are dominated by bright line emission regions (usually high density and 
high temperature regions) within each cell while the average temperature of each cell is lower. 
This interpretation is consistent with the fact that the systematic deviation becomes weaker 
for analysis using the smaller aperture size [right panel). 

Figure 6 plots the fraction of baryon mass contained in those cells where Ovil (left) 
and Ovill (right) line emissions exceed the given detection limit. Three filled patterns indicate 
the different ranges of the temperature mass-weighted over the "observed" regions. For the 
same reason as described above, especially at lower temperature regions, the line emission is 
not necessarily responsible for all the mass inside the cell. Therefore, the mass fraction shown 
here does not follow the detected mass precisely and it should be rather interpreted as the 
upper limit of the "detected" fraction. These plots again indicate that more than 90 % of 
the "detected" mass has temperature higher than 10^ K and that ~ 60 % has temperature 
with lO^K < T < lO^K. Since baryons with T > 10^ K are supposed to be already detected 
through continuum X-ray emission by conventional X-ray missions, previously unexplored and 
detectable only through the oxygen emissions will be about 20 % for the nominal detection 
limit of BIOS, 10~^^erg s~^ cm~^ sr~^. It should be also noted that the "detected fraction" 
depends on the metallicity assumed here (eq. [1]), and that if higher (lower) metallicity models 
are assumed, we will have higher (lower) "detected" fraction for a given detection limit. 

For the detection of WHIM through its emission, outskirts of rich galaxy clusters are the 
most promising regions for the targeted survey. We estimate the probability that Ovill (653 
eV) and Ovil (574 eV) line emissions are brighter than a given limiting flux within a 0.5° x 0.5° 
region around rich galaxy clusters. Coma, Perseus, Hydra, A3627 and Centaurus clusters. 
Incidentally figure 3 indicates that the brightest Ovil and Ovill emitter in this simulation is 
Virgo cluster. However, considering the 2eV energy resolution of XSA on-board DIOS, the 
minimum redshift of an extra-galactic oxygen line emitter where it can be discriminated from 
the Galactic oxygen emission is Zmin — (2eV/600eV) = 0.0033. Therefore, oxygen emission lines 
of Virgo cluster at redshift z = 0.0038 will be seriously contaminated by the strong Galactic 
oxygen emission (McCammon et al. 2002). This is why we do not attempt the mock observation 
of Virgo in the present analysis. 

Figure 7 shows such probabilities for those regions whose separation r^ from cluster cen- 
ters is Ts < l/i~^Mpc, l/i"^Mpc < Ts < 2/i"^Mpc, and 2/i"^Mpc < rg < 4/i"^Mpc. Near the central 
regions of galaxy clusters, Ovil emission cannot be detected so frequently as Ovill emission 
because the emissivity of Ovil rapidly drops at high temperature T > lO^K. Approximately 
(20-30)% area of the outskirts of known galaxy clusters {lh~^Mpc < rg < 4/i~^Mpc) in the local 
universe (Coma, Hydra, Centaurus, A3627, and Perseus) exceeds the nominal detection limit 
of DIOS (10~^^erg s^^ cm~^ sr~^) for OviI and OviII emissions. 
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4. Examples of Mock Spectra with BIOS 

Let us provide more specific examples of mock observation of the local universe with 
Dies' including the Coma region (Figs. 8 and 9), the filament in front of A3627 (Figs. 10 and 
11), and a clump roughly along the line of sight toward (simulated) Coma (Figs. 12 and 13). 
They are intended to represent rich clusters, a filament extending from clusters, and a relatively 
small clump of WHIM, respectively. 

4.1. Simulated Coma 

Figure 8 depicts four maps toward the simulated Coma cluster at redshift z = 0.023; 
X-ray flux in 0.5-2 keV band, temperature weighted by emission in the 0.5-2.0 keV band, and 
Ovil (574 eV) and Ovill (653 eV) line intensities. The corresponding mock BIOS spectra for 
the selected 3 regions (labeled 1, 2, and 3 in Fig 8) are shown in Figure 9 (500eV <E < 700eV). 
The field-of-view of those regions is chosen as 1° x 1°. 

As in Paper 1, we consider the CXB and the Galactic line emission as contaminating 
sources, and modeled according to Kushino et al. (2002) and McCammon et al. (2002), respec- 
tively. Figure 9 represents the residual spectra after the CXB and the Galactic line emissions are 
statistically subtracted. The flux intensity at an energy range from E to E + AE is calculated 
superposition of spectra for SPH particles: 

r{E+AE){l+Zi) 

X / P{E',Ti,Z,)dE', (4) 

where the summation performed over all SPH particles in the selected angular region. The 
template spectrum, P{E,T,Z), for temperature T and metallicity Z calculated using the SPEX 
software version 1.10 (http://www.sron.nl/divisions/hea/spex). 

The exposure times Texp for each region are shown in each panel. The region 1, the 
center of Coma cluster, exhibits strong continuum and Ovill emission signiflcantly larger than 
the CXB flux, while Ovil emission is faint because of its high temperature (T ^ 10^ K). In 
contrast, the continuum emission from region 2, the outskirts of Coma cluster, is much fainter 
than the CXB, and both of Ovil triplet and Ovill emission lines are not detectable. This is 
because the temperature in this region is still high (T ^ lO^K) and the emission measure is too 
low to be observable. Finally for region 3 which contains a small clump with lower temperature, 
Ovill line emission can be detected with sufficient signal-to-noise ratio if Tcxp = lO^sec, while 
the detection of Ovil triplet lines requires longer exposure time. Therefore it is not easy 
to identify the oxygen line emissions from outskirts of rich clusters, like Coma, with BIOS 
primarily because the temperature is too high for the appropriate oxygen ionization states. 
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4-. 2. the filament in front of simulated A3627 

Consider next the filament in front of simulated A3627 as an example of prominent 
filamentary structure with moderate temperature. Figure 10 shows the maps similar to Figure 8 
but for the sky area toward the east-side of simulated A3627 cluster. Note that A3627 is located 
at SGL = 190° and SGB = 7°, outside of the plots. The filamentary structure in this plot is 
located in front of A3627 as is shown in Figure 2. Figure 11 plots mock spectra for 3 regions 
marked in Figure 10. For those three regions, the continuum spectrum of the CXB contribution 
is significantly larger than that of the WHIM. Nevertheless Ovil and Ovill emission lines from 
this filament are very prominent in the subtracted spectra. For region 1, Ovill emission shows 
up at photon energy E = 641 eV corresponding to redshift z = 0.018, as well as Ovil triplet lines 
at the same redshift. In the spectra of the regions 2 and 3, Ovill and Ovil triplet lines are seen 
at z = 0.012. All these line emissions come from substructures or their outskirts in the filament. 
Contrary to rich clusters, the gas temperature in this filamentary structure lO^K <T < 10''K is 
indeed suitable for identifying Ovil and Ovill lines. In other words, small galaxy groups and/or 
outskirts of moderate clusters, which are expected to have gas temperature lO^K < T < lO^K, 
are potentially the most important targets in searching for WHIM signatures with DIOS, which 
will unveil the large-scale filamentary structure which cannot be identified with conventional 
X-ray observations. 

4-3. a clump along the line of sight toward Coma region 

In plotting Figure 2, we noticed the presence of a small clump of oxygen emitting WHIM 
along the line of sight toward Coma region {z ~ 0.01). Figure 12 is the projection maps on 
the sky toward the small clump. The center of the clump is located at SGL = 89° and SLB = 
—3°, whereas the simulated Coma (Fig. 8) is at SGL = 90° and SGB = 7°. Thus the physical 
separation between the centers of the clump and the line-of-sight toward Coma corresponds 
to ~ 5h~^Mpc at z = 0.01. Since the initial condition of the simulation is generated from 
real galaxy distribution but smoothed over 5/i~^Mpc, such a level of misalignment may be 
reasonably accounted for by an inevitable uncertainty/fiuctuation, and it is indeed possible 
that they are located almost along the same line of sight. As discussed in detail in the next 
section, this is why we are interested in this particular clump. 

Figure 13 shows the mock spectra for the selected three regions (labeled 1, 2, and 3) 
in Figure 12. For all of the three regions, we have strong Ovill line emissions at E = 647 eV, 
which correspond to redshift z = 0.009, and the Ovil triplet emission lines are also located at 
the same redshift. The peak amplitude of the X-ray surface brightness (0.5 - 2 keV band) of 
this clump is a factor of 100 times lower than that of Coma cluster. On the other hand, Ovill 
emission intensity is almost the same as Coma cluster, and Ovil intensity is even larger. 
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5. Excess emission and absorption features around the simulated Coma Cluster 

Finoguenov, Briel, & Henry (2003) reported the XMM-Newton observation of peripherals 
of Coma cluster (~ 40 arcmin apart from the center of Coma). They detected Ovil and Ovill 
emission from a possible WHIM clump with redshift estimated as z = 0.007 ± 0.004 ± 0.015 
(the first and second errors indicate statistical and systematic ones, respectively) just in front 
of Coma cluster. Since the error of the estimated redshift is large, one cannot rule out a 
possibility that the oxygen emission is associated with Coma cluster {z ~ 0.023 or ~ 70/i~^Mpc 
in distance) or our Galaxy. Nevertheless it is quite interesting that our simulated local universe 
indeed has a possible counterpart for the WHIM at an approximately right redshift z ~ 0.01. 
Therefore it is tempting to carry out more quantitative comparison between the simulated and 
the real universes. 

Figure 14 depicts the profiles of gas overdensity (top), metallicity {middle) and mass 
weighted temperature {bottom) along the line-of-sight toward the three regions (1, 2, and 3). 
These plots suggest that the simulated clump at ~ 32/i~^Mpc in distance {z ^ 0.01) is charac- 
terized by an overdensity of 6 = 10^ — 10'^, metallicity oi Z = 0.03 — O.IZq, and temperature of 
kT = 0.1 — 1.0 keV. These physical parameters are also roughly consistent with those estimated 
by Finoguenov, Briel, & Henry (2003) (see their Table 1), though the simulated temperature 
is slightly higher than their estimation. 

In order to see if we can reproduce their observed spectra, we artificially align the centers 
of the simulated clump and the simulated Coma cluster. Then we compute the mock spectra 
of the region 4 indicated both in Figure 8 and Figure 12. The size of the region 4 is chosen as 
12' X 12' so as to match the observed region of Finoguenov, Briel, & Henry (2003). Then we 
add those spectra assuming that both of them are located along the same line-of-sight ~ 40' 
apart from the center of the simulated Coma cluster. Figure 15 plots the spectrum combining 
the two regions. This should be compared with Figure 2 of Finoguenov, Briel, & Henry (2003) 
except the fact that we do not take account of the detector response function. At photon 
energy E <1 keV, the emission from the clump exceeds that from Coma cluster and the total 
spectrum exhibits 2-3 times higher flux at that energy range. This is qualitatively consistent 
with the results by Finoguenov, Briel, & Henry (2003). 

We overlay in a long dashed curve the observed spectrum of "Coma- 11 fleld" of 
Finoguenov, Briel, & Henry (2003) on the basis of the physical parameters estimated by them. 
The resulting curve is systematically shifted toward the lower energy side because the estimated 
temperature of the observed WHIM is kT = 0.24 keV while our simulated clump has a higher 
temperature kT = 0.62 keV. Thus we also attempt to calculate the spectrum of region 4 in 
Figure 12 after artiflcially replacing the gas temperature of the simulated clump (0.67 keV) 
by the observed value (0.24 keV). The resulting spectrum of the clump is shown in Figure 13. 
Now it is quantitatively consistent with the observed spectrum within a factor of two. These 
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results lend interesting, albeit indirect, support for the interpretation that the simulated clump 
corresponds to what is reported by Finoguenov, Briel, & Henry (2003). 

Consider next a possibility to detect the clump in front of Coma cluster through its 
absorption feature. The detection of metal absorption lines using the grating spectrometers 
enables to determine the redshift of the clump more accurately. Thus it serves as a complemen- 
tary and independent verification of the presence of the WHIM clump. Fujimoto et al. (2004), 
for instance, reported the 2.7a detection of Ovill absorption features aroud the outskirts of 
Virgo cluster. Indeed XMM-Newton is scheduled to observe the Seyfert I galaxy X-Comae lo- 
cated ~ 30 arcmin from the center of Coma cluster with 300 ksec exposure (J. P. Henry, private 
communication) in order to detect the possible absorption feature due to the WHIM associated 
with Coma cluster and the clump suggested by Finoguenov, Briel, & Henry (2003). As men- 
tioned in the above, the temperature of the clump (0.24 keV) reported by Finoguenov, Briel, & 
Henry (2003) is lower than the simulated temperature (0.67 keV). So again we try to compute 
the expected significance of the Ovil and Ovill absorption lines for two cases; in the first case, 
we adopt the simulated temperature as it is, and in the second case, we assume that the gas 
clump has an isothermal temperature of 0.24 keV. 

For the first case, the Ovil and Ovill column densities for the three regions (1,2, and 3) 
amount to (1.9 — 5.7) x 10^^ cm~^ and (2.0 — 5.4) x 10^^ cm~^, respectively. Here, we calculated 
Ovil and Ovill ionization fractions with the CLOUDY code (Ferland et al. 1998) assuming that 
the baryon in the clump is in coUisional ionization equilibrium. If we assume photo-ionization 
equilibrium adopting sum of the UV background estimated by ShuU et al. (1999) and the X-ray 
background by Miyaji et al. (1998) as the background photon field, the Ovil and Ovill column 
densities for those regions are (1.4 — 5.1) x 10^^ cm~^ and (2.3 — 5.9) x 10^^ cm~^, respectively. 
Considering the X-ray flux of X-Comae of (1.6 — 3.9) x 10~^^ erg s~^ cm~^ at (0.1-2.0) keV 
band and the simulated Ovil and Ovill column densities described above, the significance of 
Ovil (574 eV) absorption line obtained by XMM-Newton observation with exposure of 3 x 10^ 
seconds is expected to be (0.75-3.5)cr. Therefore, we may marginally detect the absorption 
feature of the clump. On the other hand, the detection of Ovill (653 eV) absorption will be 
not significant (less than 2a even in the most optimistic case). We simulate the XMM-Newton 
RGS observation of X-Comae through the gas clump in the region 1 in the Figure 12. The 
resulting spectrum is shown in the upper panel of Figure 16. Here we set the exposure time to 
300 ksec. At wavelength A = 19. 2A and 21. 8A, we can see very weak absorption features due 
to Ovill (653 eV) and Ovil (567 eV) around redshift z ~ 0.01, respectively. As expected from 
the computed column density of OviI and OviII, the absorption features are too weak to be 
detected. 

For the second case, on the other hand, the resulting column densities are (0.25 — 2.6) x 
10^^ cm~^ and (0.60 — 5.8) x 10^^ cm~^ for Ovil and Ovill, respectively. The expected signifi- 
cance of Ovil and Ovill absorption lines are (0.69 — n)a and (1 — 22)cr, respectively. Therefore, 
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in this case, in contrast to the first case, we expect the firm detection of both Ovil and Ovill 
absorption hnes. The lower panel of Figure 16 shows the simulated XMM-Newton RGS spec- 
trum for the second case. All the model parameters except the gas temperature are the same 
as those adopted in the upper panel. In this case, Ovil (574 eV) and Ovill (653 eV) absorption 
lines will be clearly detected with sufficient significance. 

6. Summary and conclusions 

We have explored the possibility to locate the Warm-Hot Intergalactic Medium in the 
nearby universe with a dedicated soft X-ray mission to search for cosmic missing baryons, 
DIOS. A constrained hydrodynamic simulation of the local universe has enabled us to construct 
realistic mock spectra assuming the energy resolution of ~ 2eV and typical exposure times of 
~ 10^ sec. We find that, while Virgo cluster is too close to our Galaxy in redshift space, and 
its oxygen emission may be severely contaminated by the Galactic one, several filamentary 
structures around the other rich clusters in the local universe are plausible and promising 
targets for the WHIM search with DIOS. Thus their targeted survey plays a complementary 
role to the blind survey for WHIM located at 0.0 < z < 0.3 as discussed in Paper I. In those 
regions both Ovil and Ovill emissions are likely to be identified simultaneously, and one can 
use their line ratio as a diagnostics for the nature of the observed WHIM, in particular its 
temperature. We also find that we can detect detect Ovil or Ovill emission in approximately 
(20-30)% area of the outskirts of known galaxy clusters in the local universe (Coma, Hydra, 
Centaurus, A3627, and Perseus). 

Interestingly our simulated local universe indeed has a possible counterpart for the 
WHIM in front of Coma cluster as reported by Finoguenov, Briel, & Henry (2003). While 
the temperature of the simulated clump is a bit higher than the observational estimate, the 
simulated mock spectra superposing the clump and the outskirts of the simulated Coma well re- 
produce the observed feature if we set the gas temperature to the observed value. Furthermore, 
we have explored a possibility to detect the clump through its absorption feature in the spec- 
trum of a bright Seyfert I galaxy, X-Comae, behind Coma cluster, and we may detect the Ovil 
absorption line due to the clump with the scheduled XMM-Newton observation of X-Comae. 

We thank Yehuda Hoffman for calling our attention to their earlier paper (Kravtsov, 
Klypin & Hoffman 2002) addressing the oxygen emission signatures from the constrained local 
universe simulations, and also for useful discussion. Numerical computations presented in this 
paper were carried out at ADAC (the Astronomical Data Analysis Center) of the National 
Astronomical Observatory, Japan (project ID: mky05a). KY and TF acknowledge support 
from the Japan Society for the Promotion of Science. This research was supported in part 
by the Grants-in-Aid by Japan Society for the Promotion of Science (14102004, 14204017, 
15340088, 15740157, 16340053). 
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Table 1. Simulated and observed properties of nearby clusters. 



Name 


dKjL [aegj 


otj-t) [aegj 


redshift 


l\/r.^oo ri ol4 n/f 1 

iviass [iU Mq\ 


Temperature [keV] 
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Q QO 

o.oz 


O OOQ 

U.Uzo 


A OQ* ^ \ 

4.y» (r < U200) 


Q ot ( ^ ^ 90'\ 

o.U' (^r < zu j 


L-oma (^sim.j 


yu.z 


7 1 O 
( .iU 


O OO/l 

U.UZ4 
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/'^ ^ 90M 
D.z (^r < zU j 
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Q 9^ 
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1 n7 
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11/1 
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7 0/1 
( .U4 


O 01 
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4.U'' (^r < 4U j 
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D.o" (^r < zU j 


AoDZ ( (^sim.j 


1 QO 

ifsy 


7 QO 


O 01 7 


1. 1 < 4U j 


Q Q ^ 90'\ 

0.0 (^r < zU j 


Hydra (obs.) 


1 

ioy 


Q7 


O 01 Q 
U.Uio 


Z.8 (,r<K2ooj 


Q Q# ^ oo'\ 
0.0^ (^r < zu j 


Hydra (sim.) 
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-36.0 


0.011 


2.6 (r < i?9nn) 


3.4 (r < 20') 


Centaurus (obs.) 


156 


-11.4 


0.011 


1.0** (r < 50') 


4.0# [r < 10') 


Centaurus (sim.) 


159 


-4.5 


0.013 


2.4 (r < 50') 


4.2 (r < 10') 


Perseus (obs.) 


348 


-14.1 


0.018 


9.1* (r<i?20o) 


6.7II (r < 20') 


Perseus (sim.) 


348 


-11.8 


0.016 


8.74 (r < i?20o) 


6.9 (r < 20') 



* Girardi et al. (1998) 
t Arnaud et al. (2001) 

* Shibata ct al. (2001). 
§ Tamura ct al. (1998) 



II Dc Grandi & Molendi (2002) 

# Purusho ct al. (2001) 

** Ettori, Dc Grandi, & Mclondi (2002) 
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Table 2. DIOS Mission summary 



Four-reflection X-ray telescope 
FOVxcfl:cctive area 
angular resolution 
X-ray imaging spectrograph 
energy range 

energy spectral resolution 
size of detector 
number of pixels 
FOV 

Satellite system 



Sn =100 cmMeg^ at 0.6 keV 
3 arcmin 



O.SkeV <E< IkeV 
2eV 

> 10mm X 10mm 

-16x16 

- 1° X 1° 



orbit lifetime > 1 year 

position control accuracy < 0.5 arcmin 

weight < 400 kg 



Table 3. Detection limit of line emission for the DIOS mission at photon energy E = 600 eV. 

Exposure [s] limiting flux for S/N = 10 limiting flux for S/N = 5 
erg s^^ cm^^ sr^^ erg s^^ cm^^ sr^^ 

10-* 3.7x10-1" 1.2x10-1" 

10^ 6.0x10-" 3.0x10-" 

10^ 1.4x10-" 6.6x10-12 




Fig. 1. All-sky map of hydrogen column density in the simulated local universe (in the supergalactic 
coordinate). Counterparts of several clusters in the real universe are labeled in the map. 
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figure2.jpg 



Fig. 2. Projected maps of gas density (upper left), temperature [upper right), and metallicity (lower 
left) in the simulated local universe on the supergalactic plane. The thickness of the slice is 9ft.^^Mpc 
(— 4.5/i"^Mpc < SGZ < 4.5/i~^Mpc). The distance is measured in real space (not in redshift space). 
Labels "H-C" and "P-P" indicate the Hydra-Centaurus supercluster and the Pisces-Perseus supercluster, 
respectively. 
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figure3a.jpg 



figure3b.jpg 
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figure3c.jpg 



figure4a.jpg 



figure4b.jpg 



Fig. 4. Projected surface brightness maps of Ovii (upper) and Ovili (lower) emission lines in 
the simulated local universe on the supergalactic plane. The thickness of the slice is 9/i^^Mpc 
(-4.5/i~^Mpc < SGZ < 4.5/i~^Mpc). The distance is measured in redshift space (0 < z < 0.03). 
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Fig. 5. Relation between the ratio of Ovii (574 eV) and Oviii (653 eV) line intensity and baryon tem- 
perature for the regions where both the surface brightness of Ovii (574 eV) and Oviii (653 eV) exceed the 
threshold 10~^^ and 10~^^ erg s^^ cm^^ sr^^ simultaneously. The three contours from inside to outside 
indicate that the 30, 60, 90% of the all regions arc enclosed inside them. The size of observed area is set to 
ASGL = ASGB = 0.5° {left) and 0.25° {right), and the redshift depth is set to Az = 0.0033 in both panels. 
The solid line indicates the theoretical relation under the assumption of coUisional ionization equilibrium. 
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Fig. 6. Baryon mass fraction detected through Ovii (574 eV) (right) and Oviii (653 eV) (left) emission 
hnes. Contributions of baryons with T < 10^ K, lO^K <T < lO^K, and T > lO'^K are shown separately. 
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Fig. 7. Cumulative distribution of the oxygen flux for a 0.5° x 0.5° region separated from rich galaxy 
clusters; {Left) for Oviii, and (Right) for Ovii. Different lines correspond to the different regions of radius 
Ts away from the center of those clusters; < l/i~^Mpc (solid), lh~^ < rg < 2/i~^Mpc (dotted), and 
2h^^ <rs < Ah^^Mpc (dashed). 
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figure8.jpg 



Fig. 8. Maps of soft X-ray (0.5-2 kcV) {upper left), emission weighted temperature {upper right), Ovii 
(574 eV) {lower left) and Oviii (653 eV) {lower right) toward simulated Coma cluster. 
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Fig. 9. Spectra of the 3 regions marked in Figure 8 after the CXB and the Galactic emission arc sub- 
tracted. 
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figurel0.jpg 



Fig. 10. Maps of soft X-ray (0.5-2 kcV) {upper left), emission weighted temperature {upper right), Ovii 
(574 eV) {lower left) and Oviii (653 eV) {lower right) toward tlie simulated A3627 cluster. 
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Fig. 11. Spectra of the 3 regions marked in Figure 10 after the CXB and the Galactic emission are 
subtracted. 
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Fig. 12. Maps of soft X-ray (0.5-2 kcV) {upper left), emission weighted temperature {upper right), Ovii 
(574 eV) {lower left) and Oviii (653 eV) {lower right) toward tlie elump in front of Coma cluster. 
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Fig. 13. Spectra of the 3 regions marked in Figure 12 after the CXB and the Galactic emission are 
subtracted. 
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distance [h~'Mpc] 



Fig. 14. Profiles of ovcrdcnsity (top), metallicity (middle), and mass-weighted temperature (bottom) for 
the 3 regions (labeled 1, 2, and 3) marked in Figure 12. 
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Fig. 15. Simulated X-ray spectra for tiie region 4 in the outskirts of simulated Coma cluster (Figure 8) 
and in the simulated clump in front of Coma (Figure 12). The spectrum using the physical parameters of 
the filament reported by Finoguenov, Briel, & Henry (2003) is also shown for comparison in the long-dash 
line. The dotted-dash line indicates the spectrum for region 4 in Figure 12 but the baryon temperature is 
set to fcr= 0.24 keV. 
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Fig. 16. Simulated XMM-Newton RGS spectra of absorption lines in the spectrum of X-Comae induced 
by the simulated gas clump in front of Coma cluster along the line-of-sight toward region 1 in Figure 12. 
In the upper panel, the simulated temperature of the clump is adopted as it is, while we assume the 
temperature of 0.24 keV estimated by Finoguenov, Briel, & Henry (2003) in the lower panel. Note that 
the absorption-like features at 20. SA in both panels are not physical but simply due to the gap of CCDs 
on the RGS detector. Solid lines are the theoretically expected spectra and points with error-bars indicate 
the simulated photon counts. 
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